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Разработана методика расчетного моделирования нейтронно-физического и теплогидравлического состояния канала водоох- 
лаждаемого реактора с графитовым замедлителем в быстропротекающих процессах. Произведен расчет гипотетической ава- 
рийной ситуации, возникающей при отказе приводов насосов и падении давления в контуре многократной принудительной 
циркуляции, которое вызвано частичным разрывом трубопровода. 


Введение 

Технические системы большой сложности и по- 
вышенной мощности неизбежно создают определен- 
ную степень риска для человека и окружающей среды. 
Можно лишь сократить до минимума вероятность тя- 
желой аварии и стремиться уменьшить ее возможные 
последствия. Проектом АЭС рассматриваются воз- 
можные аварии, для которых известны гипотетиче- 
ские исходные события и конечное состояние. Для 
них создаются системы безопасности, ограничива- 
ющие последствия установленными пределами. 

Расчетные пакеты прикладных программ явля- 
ются основным инструментом в анализе тепловых и 
гидродинамических процессов канальных реакто- 
ров, в том числе при частичных разрывах контура 
многократной принудительной циркуляции РБМК 
(реактор большой мощности канальный). Характер- 
ной чертой исследований безопасности РБМК при 
авариях с потерей теплоносителя является широкое 
использование зарубежных пакетов прикладных 
программ К Б БАР АТНБЕТ и др. [1]. Они созданы 
применительно к анализу безопасности реакторов 
типа Р\ѴК (ргешгеб тіег геасіог), поэтому анализ 
аварийных режимов РБМК на их основе имеет 
свою специфику. Прежде всего, встает вопрос о воз- 


можности применения кодов в условиях, опреде- 
ленных особенностями канального реактора с гра- 
фитовым замедлителем [2]. В кодах РЕБАР и АТН- 
БЕТ в явном виде не моделируется динамика паро- 
образования и динамика двухфазного потока. Для 
Р\ѴК и ВВЭР ( водо-водяной энергетический реактор) 
в этом нет необходимости. Поэтому при использо- 
вании указанных пакетов прикладных программ 
нет возможности определить, какой эффект прева- 
лирует - отрицательная обратная связь за счет Доп- 
плер-эффекта или положительная обратная связь, 
обусловленная парообразованием. 

Цель работы и методика решения задачи 

Целью настоящей работы являлась разработка 
методики расчетного моделирования нейтронно- 
физического и теплогидравлического состояний 
канала реактора в быстропротекающих процессах. 

Быстропротекающие процессы в реакторах мо- 
гут инициироваться нештатными ситуациями, в 
частности, неконтролируемыми изменениями: 

- расхода и давления теплоносителя в контуре 
многократной принудительной циркуляции; 

- реактивности и, соответственно, нейтронной 
мощности. 
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Основным принципом, используемым в мето- 
дике, является явный учет обратных связей между 
возмущением, вносимым в систему, и откликом 
системы на возмущение. В качестве возмущения 
выступает изменение реактивности, в качестве от- 
клика - нейтронная мощность. Она определяет те- 
плофизическое состояние активной зоны, измене- 
ние которого, в свою очередь, обусловливает сте- 
пень отклонения реактора от критического состоя- 
ния, т.е. реактивность. Такая модель системы отно- 
сится к классу самосогласованных- замкнутых. 

Связь между реактивностью и нейтронной 
мощностью описывается в терминах функции ли- 
нейного отклика [3]. 


Быстропротекающий процесс дробится во време- 
ни; в пределах каждого интервала времени путем ор- 
ганизации итерационного процесса решается самос- 
огласованная нейтронно- и теплофизическая задача. 

На рис. 1 представлена структура модели систе- 
мы и алгоритм вычислений, которые составляют 
основу пакета прикладных программ. 

Это частный случай, поскольку обратная связь в 
нем образована с учетом изменения физических 
свойств теплоносителя - воды (пароводяной смеси). 

На первом временном интервале А г, начиная с 
4=0, для фиксированного стационарного уровня 
мощности и 0 =1 отн. ед., при заданном расходе тепло- 



Рис. 1 . Структура модели и организация вычислений 
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носителя С 0 и при неизменном давлении Р 0 прово- 
дится расчет теплофизического (ТФ) и гидрогазоди- 
намического (ГГД) состояния технологического ка- 
нала (ТК). Рассчитываются изменение температуры 
и скорости теплоносителя при движении в ТК, ко- 
эффициенты теплоотдачи с поверхности твэла, про- 
водится оценка критических тепловых потоков. 

В результате определяется температурное со- 
стояние элементов ТК по аксиальной координате 
I, в т.ч. плотность потока теплоносителя (ТН) - 
у Н;0 (ф. Затем проводится расчет параметров диа- 
батного потока, результатом которого является 
значение истинного объемного паросодержания, 
и, следовательно: 

- площади проходного сечения ТК, занятого па- 
ровой О 5н г0 и жидкой °У Ні0 фазами; 

- плотности молекул Н,0 в паровой °ун 2 о и жид- 
кой Уцр фазах. 

Далее следует задание новых значений расхода 
и давления - С ] и Р х . Новое значение нейтронной 
мощности полагается неизменным п\=щ. Цикл 
расчетов (ТФ и ГГД)— > расчет параметров диабат- 
ного потока повторяется. В результате (/=1) опре- 
деляются новые значения ‘.5” Н20 , ‘5^0, Ѵн 2 о и Ѵн 2 о- 
Приращения ДУнр, АіУн,о и изменения значе- 
ний Д'нр, 5 н 2 о, определяют значение приращения 
реактивности Ар,. 

Затем определяется «истинное» значение ней- 
тронной мощности й[ в формализме, использую- 
щем понятие функции линейного отклика. 

Далее следует изменение расхода и давления: 0 2 
и Р 2 . Новое значение нейтронной мощности пола- 
гается неизменным л,* =й і- Цикл расчетов [(ТФ и 
ГГД)— > расчет параметров диабатного потока] пов- 
торяется. В результате (/=2) определяются новые 
значения 2 5” Н2 о, 2 ^н 2 о, Ун 2 о и 2 уЁ і0 и следует нейтрон- 
но-физический расчет. 

Начиная с /= 1 последовательность вычислений 
остаётся неизменной и повторяется. Результатом 
вычислений является динамика нейтронно-физи- 
ческих и теплогидравлических параметров ТК. 

Последовательность вычислений 

при определении параметров диабатного потока 

Расчет по соотношениям, которые используются 
при определении параметров диабатного потока [4], 
проводится в последовательности, изложенной ниже. 

В соотношениях приняты следующие обозначе- 
ния: р - давление, МПа; 1 - температура, °С; С - мас- 
совый расход, кг/с; Р— проходное сечение, м 2 ; сі- ди- 
аметр твэла, мм; - плотность теплового потока, 
МВт/м 2 ; ( рсо)=6/Р - массовая скорость, кг/м 2 -с. 

Из подпрограмм-таблиц по заданным р и / 
определяются: г - удельная теплота парообразова- 
ния (теплота испарения), Дж/кг; И - удельная эн- 
тальпия потока, Дж/кг; /Г - для жидкой фазы, 
Дж/кг; к для паровой фазы, Дж/кг; ѵ - удельный 
объём потока, м’/кг; ѵ'— для жидкой фазы, м’/кг; 


р — 1/ѵ'— плотность жидкой фазы, кг/м 3 ; ѵ"— для 
паровой фазы, м 3 /кг; р'—І/ѵ"— плотность паровой 
фазы, кг/м 3 ; ѵ' - коэффициент кинематической 
вязкости для жидкой фазы, м 2 /с; для заданной ( 
определяется р кр - критическое давление, МПа. За- 
тем проводятся вычисления: 


1) со 0 =рсо/р'- скорость циркуляции; 

2) я=р/р ѵ ; 

3) коэффициент скольжения: 

5=1+2,27(1 -тг 2 )(р/р®) 0 ' 7 ; 

4) средняя энтальпия потока в точке начала ин- 
тенсивного роста истинного объёмного паросо- 
держания: 


рсо 


гр”ѵ' 



5) АА = (А + А'-2А 0 )ехр^-2-^— ^- 1 ; 

6) массовое расходное паросодержание: 

х = (А- А' + АА)/(г + А/?); 

7) граничное значение истинного объёмного па- 
росодержания: 

ср р = ( 1+0, ЗЗЗ.?) -1 ; 

8) граничное значение массового расходного па- 
росодержания: 

х=(1+0,ЗЗЗрУрУ; 

9) <р 0 =1,17<У// 15 (рй>)° 15 ; 

10) х 0 =-О,537 ф 1 (р/ра>)'»; 

Истинное объёмное паросодержание потока (р: 

1 1) ^=0 при х<х 0 (п. 6, п. 10), 

(р = ср 0 (1-х/х 0 ) 1,35 при х 0 < х < 0, 

<Р = ( Р 0 +—(<Р р -<Р 0 ) при 0 < х< х р , 
х р 

<Р = <Р Ѵ при х>х р . 

Площадь проходного сечения ТК, занятого па- 
ровой фазой - 5н р=(рР, занятого жидкой фазой - 
5 Н! о=Т-5н !0 . 


Расчет отклика нейтронной мощности 
на возмущение реактивности 

Выражения для концентрации ядер-предше- 


ственников запаздывающих нейтронов і- 

ои груп- 

пы при: 




г=0 

'■° р я,. ' 
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— сі п Ал а (1 
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0<т^Дг 
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Здесь а,=Д/Д - относительный выход запазды- 
вающих нейтронов по группам; Ап=п ] —п {] - изме- 
нение мощности за время А г. 


Функция линейного отклика 


сх 


Подстановка выражений ( 1-3) в Р = 1 - X позволя- 
ет получить следующее выражение для изменения реак- 
тивности при заданном законе изменения мощности: 

- Ап тт / ч 

Р= — Л»- 
п(т) 

При указанной подстановке учтено, что 

І>,= 1. 


Функция (7 Дг (т) называется функцией линейно- 
го отклика [3]. 

Эта функция задается следующими соотно- 
шениями при: 


0<т<Дг 


^)=х 


а / (1-е~ я,г ) 
Я,. А г 


і~\ я,.дг 

Из (3) получаем следующее значение конечного 
уровня мощности реактора л,, которое достигается 
после ввода реактивности р за время Ат 


1± — к 

^дДАт) 


В нашем случае возмущение реактивности об- 
условлено изменением значения истинного объё- 
много паросодержания (изменением концентра- 
ции молекул воды) 


Ар = А и н 


сір 

Ѵ^Ѵо у 


Для реактора с отражателем необходимо учиты- 
вать «эффективные добавки» за счёт отражателя. В 
этом случае, при расчёте геометрического параметра, 
фактические полуразмеры зоны увеличиваются на ве- 
личину эффективной добавки 8. Н-Н+28', Я-К+8. 

Для уран-графитовых реакторов с графитовым 
отражателем 

8 = 1,2Ь отр іЪ-^, 

^ отр 


где Ь отр - длина диффузии в отражателе (графите), 
Ь ощ =5 7 см, Т - толщина отражателя, см. 

Для реактора РБМК-1000 толщина отражате- 
лей: бокового - 75 см; верхнего и нижнего - 50 см. 


где г с =312,5 см 2 - возраст нейтронов в 

графите, 8 Ш , 8 Ш - площади поперечного сечения 
ячейки и ТК, 8 КМ - конструкционных элементов ТК, 
А н ,о _ проходного сечения под пароводяную смесь. 

Вероятность избежать резонансного захвата 
определяется соотношением: 


ср = ехр 


(к т К и у[Кпе і + 0,73- пВ^е^) 


где и <рі_ 5 ф - замедляющая способность замедли- 
теля и «фиктивного блока»; Д 3 и 8 Ф - площадь замед- 
лителя и фиктивного блока; К ѵ - радиус уранового 
блока; к т - температурный коэффициент; п - число 
стержней в пучке; Я - радиус пучка; е, - пористость 
блока по урану-238. Коэффициент к т имеет вид: 

к т =0,775(1 + 17, 5 -10 -3 ^7), 


где К и - средняя температура урана, К. 

Изменение концентрации молекул воды наибо- 
лее сильно влияет на коэффициент использования 
тепловых нейтронов Ѳ , поэтому при дифференци- 
ровании формулы четырех сомножителей считаем 
остальные сомножители константами: 


Р = 1- 


1 + Ь 2 В 2 


т)8(рѲ ■ ехр(— 2?~г ) 


Р = 


1 


ѵ реср -ехр(-В 2 т ) у 

' 1 


А + Ь 2 В 2 ) _ 
' Ѳ ’ 


?!8(р -ехір(- В т) 


•(-^(і + Т-в 2 )), 

V 


сір __ СІѲ \ + Ь 2 В 2 

4,0 ^ ,7 н 2 0 Ѳ 2 г]8(р-ехр(-В 2 т) 


(4) 


Дифференцируем выражение для определения 
коэффициента использования тепловых нейтронов: 

у топя у топл 

Ѳ = ^~ 






НлО 


сІѲ 

СІП и л 


=хі 


.ѵ~ ,2г х -1 N6 х~ ' с ѵ -4 Н,0 X - ' топл \ 2 

Ѣа+Х'+Х+Ха ) 


При расчете по соотношению (4) значения Ѳ и 
определяются для 23° которое вычисляется 

по формуле: 

піП П I С' ' 

, И .О , н ,0 У н 2 о н,о7н,о 

) — \^а > Л > 

*^Н 2 0 ’ 7 н 2 о 

где (Е® ! °)* - макроскопическое сечение поглоще- 
ния в пароводяной смеси (в воде) на предыдущем 
шаге интегрирования по времени. Для воды 


= ^--0,0221, см -1 , где — плотность воды в 

У Н 2 0 

нормальных условиях, у Ні0 - при текущих значе- 
ниях температуры и давления. 


67 




Известия Томского политехнического университета. 2005. Т. 308. № 3 


Результаты расчетов и их обсуждение 

В расчетах моделировалась гипотетическая ава- 
рийная ситуация, возникающая при отказе приво- 
дов насосов и падении давления в контуре много- 
кратной принудительной циркуляции, которое вы- 
звано частичным разрывом трубопровода. 

Динамика теплового состояния технологиче- 
ского канала приведена на рис. 2. 

Расход теплоносителя за 1,8 с уменьшается в 
1,3 раза, и в этот момент возникают условия, 
необходимые для инициации местного кипения на 
наиболее энергонапряженном участке ТК. Темпе- 
ратура поверхности твэла достигает 350 °С, превы- 
шая температуру насыщения на 65 °С. За следую- 
щие 0,5 с возникают условия, необходимые для 
объемного кипения. В случае, если паросодержа- 
ние достигает граничного значения, наступает кри- 
зис теплообмена 2-ого рода: жидкая пленка высы- 
хает, вся выделяющаяся энергия идет на разогрев 
твэла, в том числе - его оболочки. Это иллюстри- 
руется ходом зависимости температуры поверхно- 
сти твэла от времени (рис. 2, зависимость 5). Тем- 
пература теплоносителя достигает температуры на- 
сыщения и далее не изменяется. Нейтронная мощ- 


ность стремительно нарастает, и через 2,5 с после 
начала падения расхода теплоносителя плотность 
теплового потока достигает критического значе- 
ния. Начиная с этого момента времени очень веро- 
ятно стремительное нарастание температуры обо- 
лочки твэла и, как следствие, ее разрушение. 

С учетом срабатывания аварийной защиты ди- 
намика процесса следующая. Время возрастания 
мощности составляет 1,5. ..2,0 с. Нейтронная мощ- 
ность сначала, вслед за ростом паросодержания, 
возрастает до 1 15 % номинального значения. После 
срабатывания защиты она быстро снижается до 
уровня остаточного тепловыделения. 

Результаты расчетов [5] показывают, что с уче- 
том того, что около 20 % запасенного в твэле тепла 
расходуется на прогрев оболочки, ее температура 
при адиабатическом выравнивании профиля тем- 
пературы может превысить 1 100 °С. Это значение 
характеризует температуру, которая может быть до- 
стигнута в авариях рассматриваемого типа. 

Оценка неопределенностей результатов, полу- 
ченных в расчетах, позволяет определить ключевые 
процессы и явления в канале водоохлаждаемого 
реактора. 


ю,м/с ц Ф , МВт/М 2 ^5, МВт/м 2 І п , : С 1, °С 



Рис. 2. Динамика теплового состояния технологического канала: 1) начало местного кипения; 2) начало объёмного кипения; 
3) I ~ температура потока; 4) Г, - температура насыщения; 5) возможный рост 1„; 6) і„ ~ температура поверхности твэ- 
ла; 7) со - скорость потока; 8) д, - плотность теплового потока; 9) возможность начала кризиса теплообмена; 
10) д кр - критическая тепловая нагрузка 
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Технические науки 


Ими следует считать: 

- запасенное тепло в ядерном топливе; 

- термодинамическое состояние системы и, в 
первую очередь - особенности течения двух- 
фазного потока, вызывающие колебания расхо- 
да теплоносителя и условий теплообмена на по- 
верхности твэла. 

Определяющими параметрами при моделиро- 
вании процесса отвода тепла, запасенного в топли- 
ве, являются: 

- термическое сопротивление зазора между то- 
пливом и оболочкой; 

- коэффициент неравномерности энерговыделе- 
ния; 
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